Experimental observation of the Anderson transition with atomic matter waves 
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We realize experimentally an atom-optics quantum chaotic system, the quasiperiodic kicked rotor, 
which is equivalent to a 3D disordered system, that allow us to demonstrate the Anderson metal- 
insulator transition. Sensitive measurements of the atomic wavefunction dynamics and the use 
of finite-size scaling techniques make it possible to extract both the critical parameters and the 
critical exponent of the transition, which is in good agreement with the value obtained in numerical 
simulations of the 3D Anderson model. 



PACS numbers: 03.75.-b, 72.15.Rn, 05.45. Mt, 



:.70.Tg 



The metal-insulator Anderson transition plays a cen- 
tral role in the study of quantum disordered systems. An 
insulator is associated with localized states of the system, 
while a metal generally displays diffusive transport asso- 
ciated with dclocalizcd states. The Anderson model [l| 
describes such a metal-insulator transition, due to quan- 
tum interference effects driven by the amount of disorder 
in the system. Starting from the "tight-binding" descrip- 
tion of an electron in a crystal lattice, Anderson postu- 
lated in 1958 that the dominant effect of impurities in 
the lattice is to randomize the diagonal, on-site, term of 
the Hamiltonian, and showed that this generally leads to 
a localization of the wavefunction, in sharp contrast with 
the Bloch-wave solution for a perfect crystal. This model 
has progressively been extended from its original solid- 
state physics scope P, d, H, 0] to a whole class of systems 
in which waves propagate in a disordered medium, as for 
example quantum-chaotic systems 0, Q and electromag- 
netic radiation 0, However mathematically simple, 
the model predicts a wealth of interesting phenomena. 
In ID, the wavefunction is always localized as recently 
observed in experiments using atomic matter waves in 
a disordered optical potential [l^, [llj ; in 3D it predicts 
a phase transition between a localized (insulator) and a 
delocalized (metal) phase at a well defined mobility edge, 
the density of impurities or the energy being the control 
parameter. 

Despite the wide interest on the Anderson transition, 
few experimental results are available. In a crystal, it 
is very difficult to obtain the conditions for a clean ob- 
servation of the Anderson localization. Firstly, one has 
no direct access to the electronic wavefunction and must 
rely on modifications of bulk properties like conductiv- 
ity 0, • Secondly, it is difficult to reduce decoherence 
sources to a low enough level. We thus engineered a 
matter-wave system that is described by an Anderson- 
like model, which allows us to probe the physics of dis- 
ordered systems in much better conditions than in con- 



densed matter physics [13[, namely: Almost no interac- 
tion between particles, weak absorption in the medium, 
no coupling with a thermal reservoir which could destroy 
localization and possibility of measuring the final quan- 
tum state of the system after a given interaction time. 
The system, the quasi-periodic kicked rotor 0, 0, Q, HE] j 
consists of cold cesium atoms exposed to a pulsed, off- 
resonant, laser standing wave. The dynamics is thus 
effectively one-dimensional along the axis of the laser 
beam, as transverse directions are uncoupled. The atoms 
interact periodically with the spatially sinusoidal poten- 
tial whose amplitude is modulated at (incommensurable) 
frequencies UJ2 and 0^3, and the corresponding Hamilto- 
nian is: 
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K cos x[\+ e cos (w2i) cos (^3^)] ^ 5{t — n) , 



n=0 



(1) 

where x is the particle position (along the laser axis), 
p is its momentum, and K is the pulse intensity. The 
pulses are short enough (compared to the center of mass 
dynamics) to be consider as instantaneous kicks. We have 
chosen normalized variables such that x is measured in 
units of the spatial period of the potential divided by 27r, 
the particle's mass is unity and time is measured in units 
of pulse period Ti . 

By taking e = in Eq. ([T]) one obtains the standard 
(strictly time-periodic) kicked rotor, a system known to 
display the phenomenon of dynamical localization 
which manifests itself by an exponential localization of 
the wavefunction in momentum space. Dynamical local- 
ization exists only when the classical dynamics of the 
kicked rotor is a chaotic diffusion. In the quantum 
case, diffusion is inhibited after some localization time by 
quantum interferences. Dynamical localization has been 
shown to be a direct analogue of Anderson localization 
in one dimension with the following correspondences: 
Localization takes place in real space for the Anderson 
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Figure 1; Phase diagram of the quasiperiodic kicked rotor, 
from numerical simulations. The localized (insulator) region 
is shown in blue, the diffusive (metallic) region is shown in red. 
The experimental parameters are swept along the diagonal 
dash-dotted line. The insets show the experimentally observed 
momentum distributions, localized in the insulator region and 
Gaussian in the diffusive (metallic) region. 



model and in momentum space for the kicked rotor; the 
"stochasticity parameter" K [sec Eq. ([Ij] corresponds to 
the ratio of hopping to diagonal energy in the Anderson 
model. The strictly random disorder in the Anderson 
model correspond to a pseudo-random potential in the 
quasi-periodic kicked rotor. Thus, the experimental ob- 
servation of dynamical localization in the kicked rotor 
14| is actually the first observation of Anderson ID lo- 
calization with atomic matter waves. 

In order to observe the Anderson transition, one must 
generalize the kicked rotor to obtain a system analogous 
to the 3D Anderson model. This can be done by in- 
troducing further temporal dependencies that make the 
Hamiltonian quasi-periodic in time, that is to take e 7^ 
in Eq. ^ with UJ2 and UJ3 incommensurate numbers. 
The resulting system has been shown to be substantially 
equivalent to the 3D Anderson model [l^ , the additional 
spatial coordinates in the Anderson model corresponding 
to the additional time dependencies in the quasi-periodic 
kicked rotor. The parameters controlling the dynamics 
are iU2 and lu^ and, more importantly, the kick strength K 
and the modulation amplitude s. In the (A', e) plane, the 
e = line corresponds to dynamically localized dynamics 
(the periodic kicked rotor equivalent to the ID Anderson 
model) and the Anderson transition takes place along a 
critical line in the plane e > 0, as shown in Fig. [TJ 

Our atom-optics realization of the kicked rotor has 
been described in detail elsewhere 3, 17, l3|. Basically, 
we cool cesium atoms in a standard magneto-optical trap, 
and, after a Sisyphus-molasses phase, we obtain a cloud 



of 10'' atoms at a temperature of 3.2 fiK. This prepares 
a sample of atoms in a thermal state whose momentum 
distribution is much narrower than the expected local- 
ization length. The atoms then interact with the optical 
potential generated by a horizontal standing wave. An 
acousto-optical modulator driven by an arbitrary-form 
synthesizer modulates the amplitude of the optical po- 
tential. One generates in this way 0.9^s-long pulses at 
Ti = 27.778 /US (corresponding to an effective Planck con- 
stant k = 2.89), to which is superimposed a modulation 
of the form Eq. ([1]), with u^l'i'K — -\/5 and W3/27r ~ vT3. 
The standing wave, of typical power 160 mW, is far off- 
resonant (7.3 GHz to red of the atomic transition, corre- 
sponding to 1.4 X 10"^ natural widths), in order to reduce 
spontaneous emission. Stimulated Raman transitions al- 
low sensitive measurements of the atomic momentum dis- 
tribution (l7j . 

The experimental values of the parameters are chosen 
according to the following considerations. Firstly, in or- 
der to prevent classical (K AM barrier) localization effects 
one must have K > 2 [l9|. For all the data presented in 
this paper, we numerically checked that the classical dy- 
namics is fully diffusive, which excludes that the observed 
transition could have a classical origin. Secondly, in or- 
der to confine the transition to a relatively narrow range 
of parameters one must cross the critical curve (Fig. [T|) 
"at a right angle"; we thus vary simultaneously the K 
and e along a line going from = 4, e = 0.1 in the local- 
ized region to if = 9, e = 0.8 in the diffusive region; the 
critical line is then crossed aX K = Kc ~ 6.6. Thirdly, 
short enough pulses must be used, so that they can be 
considered as delta pulses. Numerical simulations taking 
into account the finite pulse duration (0.9 /its) show that 
less than 1% of the atoms are sensitive to the duration 
of the pulses. Finally, decoherence processes must be 
kept as small as possible. We have identified two domi- 
nant decoherence sources in our experiment: i) Sponta- 
neous emission, which is not included in the Hamiltonian 
Eq. ID) and ii) the deviation of the standing wave from 
strict horizontality, which mixes quasimomentum classes 
and produces a stray momentum diffusion. The fact that 
numerical simulations of the "pure" quasi-periodic kicked 
rotor shown hereafter agree very well with the experimen- 
tal results (in particular for the position of the Anderson 
transition and for the critical exponent) proves that spu- 
rious effects are well under control. Moreover, we have 
checked that inclusion of these effects in the numerical 
calculations only leads to small changes for t < 150 kicks. 

In order to observe the Anderson transition we ap- 
ply a sequence of kicks to the atomic cloud and mea- 
sure its dynamics. In the localized regime, the evolu- 
tion of its momentum distribution is "frozen" after the 
localization time (^--^12 kicks) into an exponential curve 
exp (— IpI/pioc) (where pioc is the localization length). In 
the diffusive regime, the initial Gaussian shape is pre- 
served and the distribution gets broader as kicks are ap- 
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Figure 2: Temporal dynamics of the quasi-periodically 
kicked rotor. We measure the population no(t) of the zero- 
momentum class as a function of time and plotted the quan- 
tity Il^^(t) (proportional to {p^)(t)). Clearly, it tends to a 
constant in the localized regime (blue curve) and increases 
linearly with time in the diffusive regime (red curve) . Close 
to the critical point (purple curve) , it displays anomalous dif- 
fusion IlQ^{t) ~ t'^^^ , as expected from theoretical arguments 



plied, corresponding to a linear increase of the average 
kinetic energy. The insets in Fig. [1] show the experimen- 
tally observed momentum distributions, an exponentially 
localized distribution for K < Kc (blue curve), charac- 
teristic of dynamical localization, and a broad, Gaussian- 
shaped distribution for K > Kc (red curve), character- 
istic of the diffusive regime. Instead of measuring the 
full momentum distribution, it is sufficient, and much 
easier, to measure the population no(t) of the zero ve- 
locity class [2l|, as Hg ^(i) is proportional to {p'^){t) (by 
virtue of the constancy of the total number of atoms). 
We performed several experimental runs corresponding 
to points in the (A", e) plane. In each run the value of 
no(t) was recorded as the kicks were applied. We also 
recorded the background signal obtained by not apply- 
ing the Raman detection sequence, and the total number 
of atoms in the cold-atom cloud. These signals are used 
to correct the experimental data from background sig- 
nals and long-term drifts of the cloud population. Figure 
[5] shows the experimentally observed 11^^ (t) and clearly 
shows the transition from the localized to the diffusive 
regime, with an intermediate regime of anomalous diffu- 
sion. 

When one approaches the critical point from the insu- 
lator side, the localization length diverges, whereas, on 
the metallic side, the diffusion constant vanishes. How- 
ever, a strict divergence can be observed only in macro- 
scopically large samples; in small samples the divergence 
is smoothed. This fact plagued the numerical studies 
of the solid-state Anderson transition, as only a finite 
(small) lattice can be dealt with in a computer. In our 



system, a singular behavior would show up only for pro- 
hibitively large numbers of kicks, which are, in practice, 
limited to 150. To overcome this limitation, a technique 
named "finite-size scaling" [13, 



22, was introduced. 



"finite-size scaling" 
whose basic idea is to infer the scaling law allowing proper 
extrapolation of the measured localization length to an 
infinite sample. We will next show that our data obey 
the scaling laws predicted by renormalization theory, and 
this will allow us to extract the critical exponent associ- 
ated with the singular behavior at the transition. 

We have adapted the standard finite-size scaling ap- 
proach used in numerical studies of the Anderson transi- 
tion 0] assuming that, for finite interaction time, the 
quantity A(t) = n^^(i)i~^/^ is an arbitrary function 
/ ) 0, H [ill , with a scaling parameter which 

depends only on K. Using the results obtained for vari- 
ous values of t and K one can reconstruct both the func- 
tion / and the scaling parameter ^(K) (no assumption on 
the form of / is made). The result is shown in Fig. [3l^a) 
and (b) for numerical simulations and in Fig. E^c) and 
(d) for the experimental results. In both cases, the scal- 
ing hypothesis is justified by the fact that all points in 
Fig-EIa) and[3l^c) lie on a single curve (within the exper- 
imental error). The remarkable feature is the existence 
of two branches: the upper one corresponds to diffusive 
motion while the lower one corresponds to the localized 
regime, the critical point being at the tip joining the 
two branches. The scaling parameter £^{K) is plotted in 
Figs. [3l^b) andlSlJd): It represents the localization length 
in the localized regime and scales as the inverse of the 
diffusion constant in the diffusive regime. Clearly, it in- 
creases rapidly in the vicinity of the critical value Kc, 
on both sides of the transition, and it is found to be- 
have as ^ ~ \K — Kc\~^ when K Kc- This divergence 
at the critical point is a key property of the Anderson 
phase transition. However, phase-breaking mechanisms 
induced by decoherence processes smooth the algebraic 
divergence. We model such effects by introducing a small 
cutoff of the divergence due to the residual decoherence: 



f 



= a\K- Kr 



(2) 



The experimental observations and the numerical data 
have been fitted with this formula [dashed curves in Figs. 
W^) andEJd)]. We found Kc = 6.4 (very close to the 
value Kc = 6.6 obtained from the numerical simulation), 
and a critical exponent v = 1.4 ± 0.3 [l^, which is con- 
sistent with the numerical value within the error bars. 
The good agreement between the numerical simulations 
and the experimental results proves that spurious effects 
(such as decoherence) are under control. We emphasize 
that there are no adjustable parameters in our procedure, 
all parameters are determined using the atoms them- 
selves as probes. Once the existence of the scaling law 
established, it is more convenient to use a global analysis 
of the numerical data at various values of K (see refs. 



4 




Figure 3: Finite size scaling applied to the results of numerical simulations [(a) and (b)] and to the experimental results [(c) 
and (d)]. The rescaled quantity A = n^^(f)f~^'''^ is plotted as a function of time (from 30 to 10* kicks (simulation) and from 
30 to 150 kicks (experiment), for various values of K. Finite-size scaling of the raw data makes it possible to determine both 
the scaling function /, shown in (a) and (c), and the scaling parameter £,{K), represented in (b) and (d). Close to Kc, the 
behavior of ^ in (b) and (d) is well fitted by Eq. ((5)) (dashed lines), giving Kc = 6.6 (simulation) and Kc = 6.4 (experiment). 
The critical exponent \s v = 1.60 ± 0.05 (simulation) and v = 1.4 ± 0.3 (experiment). 



Mackinnon [1^, Slevin and Ohtsuki [23|). This, together 
with the fact that numerical simulations are not limited 
to 150 kicks and can be ran up to several thousands kicks, 
make it possible to obtain a more precise numerical value 
for the critical exponent u ~ 1.60 ± 0.05. Note, finally, 
that although the data displayed here concerns a partic- 
ular set of parameters we also verified experimentally the 
presence of the transition for other parameters. 

In conclusion, we have presented the first experimental 
evidence of the Anderson transition with atomic matter 
waves. The transition is characterized by a well defined 
critical point, a divergence of the localization length be- 
low the critical point (in the localized regime) and a van- 
ishing of the diffusion constant above the critical point 
(in the diffusive regime). We have determined the scaling 
laws and the critical exponent ly of the Anderson tran- 
sition, which is significantly larger than unity and very 
close to the value obtained in recent numerical exper- 
iments [12I, [23I on the Anderson model, enforcing the 
assumption [l5[ that the two systems are substantially 
equivalent. Whether this exponent is universal (i.e. in- 
dependent of the microscopic details) or not remains to 
be studied. A very interesting point is that our Anderson- 
equivalent quasiperiodic kicked rotor can be easily gen- 
eralized to higher dimensions simply by adding new in- 
commensurable frequency, which opens perspectives for 
fascinating studies of the dependence of the critical expo- 
nent on the dimension of the underlying Anderson model. 

The authors acknowledge D. Shepelyansky for bringing 
ref. [isl to their attention, and for fruitful discussions. 
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